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Abstract 

In this paper we provide an analytical procedure which leads to a system of (n — 2) 2 poly- 
nomial equations whose solutions give the parameterisation of the complex n x n Hadamard 
matrices. It is shown that in general the Hadamard matrices depend on a number of arbitrary 
phases and a lower bound for this number is given. The moduli equations define interesting 
geometrical objects whose study will shed light on the parameterisation of Hadamard matrices, 
as well as on some interesting geometrical varieties defined by them. 

1 Introduction 

Quantum information theory whose main source comes of a few astonishing features in the founda- 
tions of quantum mechanics is the theory of that kind of information which is carried by quantum 
systems from the preparation device to the measuring apparatus in a quantum mechanical exper- 
iment, see e.g. [2Zj- Defining new concepts like entangled states, teleportation or dense coding 
one hopes to be able to design and construct new devices, like quantum computers, which will 
be useful in solving many "unresolvable" problems by the classical methods. Recently the mathe- 
matical structure which is behind such miracle machines was better understood by establishing a 
one-to-one correspondence between quantum teleportation schemes, dense coding schemes, orthog- 
onal bases of maximally entangled vectors, bases of unitary operators and unitary depolarizers by 
showing that given any object of any one of the above types one can construct any object of each 
of these types by using a precise procedure. See Vollbrecht and Werner and Werner j^Hl for 
details. The construction procedure will be efficient to the extent that the unitary bases can be 
generated, and the construction of these bases makes explicit use of the complex Hadamard matri- 
ces and Latin squares. The aim of this paper is to provide a procedure for the parametrisation of 
the complex Hadamard matrices for an arbitrary integer n. More precisely we will obtain a set of 



(n — 2) 2 equations whose solutions will give all the complex Hadamard matrices of size n. Complex 
n-dimensional Hadamard matrices are unitary n x n matrices whose entries have modulus 1/y/n. 

The term Hadamard matrix has its root in the Hadamard's paper ^7j, where he gave the solution 
to the question of the maximum possible absolute value of the determinant of a complex nxn matrix 
whose entries are bounded by some constant, which, without loss of generality, can be taken equal 
to unity. Hadamard has shown that the maximum is attained by complex unitary matrices whose 
entries have the same modulus and he asked the question if the maximum can also be attained 
by orthogonal matrices. These last matrices have come to be known as Hadamard matrices in his 
honor, and have many applications in combinatorics, coding theory, orthogonal designs, quantum 
information theory, etc., and a good reference about the obtained results is Agaian pQ. 

However the first complex Hadamard matrices were found by Sylvester |24j . He observed that 
if ai, i = 0, \,...,n — 1 denote the solutions of the equation x n — 1 = for a prime n then the 
Vandermonde matrix 

/ 

1 

V 

is unitary and Hadamard. In the same paper Sylvester found a method to obtain a Hadamard 
matrix of size mn if one knows two Hadamard matrices of order m and respectively n by taking 
their Kronecker product. Soon after the publication of the paper by Hadamard the interest was 
mainly on the real Hadamard matrices such that the Sylvester contribution fell into oblivion and 
the complex Hadamard matrices have been much later reinvented in a particular case: only those 
matrices whose entries are db 1 , ± i where i = 

Nevertheless a few other problems apparently unrelated to complex Hadamard matrices were 
those connected with bounds on polynomial coefficients when the indeterminate runs on the unit 
circle. They are better expressed in terms of the discrete Fourier transform. For any finite sequence 
x = (xq,x\, . . . ,x n -i) of n complex numbers, its (discrete) Fourier transform is defined by 

n-l 

%- = ^" 1/2 E Xke 2i * kj/n j = 0,l,...,re-l 

k=0 

If the components Xk,Vk are such that \xt\ = \Vk\ = 1 for k = 0, 1, . . . , n — 1 the sequence x is called 
bi-unimodular. The existence of a bi-unimodular sequence of side re is equivalent to the existence 
of a complex circulant Hadamard matrix of side re; a circulant matrix is obtained by circulating its 
first row, in our case the components of the vector x / y/n. Now the Gauss sequence 

( e 2nr(ak 2 +bk)/n^ a> g z, a coprimeto re, A; = 0, 1, . . . , n — 1 for re odd 
1 e fc2 l7r / n ) fc = 0, 1, . . . ,n — 1 for re even 

is a bi-unimodular sequence [S]. The problem of the complete determination of all bi-unimodular 
sequences is still open, despite the problem is simpler than the parameterisation of arbitrary com- 
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plex Hadamard matrices. However this approach gave the first non-trivial examples of complex 
Hadamard matrices for n > 6. 

A step towards its solution was the reduction of the bi-unimodular problem to the problem of 
finding all cyclic n-roots [S], that are given by the following system of equations over C 

zo + zi H h z n _i = 0, 

z G zi + z x z 2 H V z n ~iZo = 0, 

< Z Z 1 Z 2 + Z 1 Z 2 Z 3 H 1" Zn-lZQZl = 0, (1) 



Z Zi ■ ■ ■ Z n -l = 1 

Note that the sums are cyclic and contain just n terms and are not the elementary symmetric 
functions for n > 4. The relation between x and z is zj = Xj+i/xj. All cyclic n-roots have 
been found for 2 < n < 8; see Bjorck and Froberg |3[7j- The formalism we will develop in the 
paper is more general showing that the parameterisation of complex Hadamard matrices is more 
complicated than the finding of all cyclic n-roots of the sytem Using our approach we find, 
e.g. when n = 6, the following matrix which is not contained in the above solutions 
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matrix that depends on an arbitrary phase. 

The parameterisation of complex Hadamard matrices is a special case of a more general problem: 
that of reconstructing the phases of a unitary matrix from the knowledge of the moduli of its entries, 
problem which was a fashionable one at the end of eighties of the last century in the high energy 
physics community, see Auberson 3 , Bjorken and Dunietz jS], Branco and Lavoura JHJ) Auberson 
et al. (I]. An existence theorem as well as an estimation for the number of solutions was obtained 
by us [121 • T ne particle physicists abandoned the problem when they realised that for n > 4 
there exists a continuum of solutions, i.e. solutions depending on arbitrary phases, result that 
was considered uninteresting from the physical point of view. In our opinion, the reason was the 
difficulty of the problem; since the experiments provide only the squares of the moduli, the first 
problem is to decide if from the experimental results, which in the best case generate a doubly 
stochastic matrix, one can reconstruct a unitary matrix, or a unistochastic matrix. Only for n = 3 
there exists a unambigous procedure. For n > 4 there are no known necessary and sufficient 
conditions to separate the unistochastic matrices from the doubly stochastic ones |29j . 

Almost in the same time the complex Hadamard matrices came out in the construction of 
some *-subalgebras in finite von Neumann algebras, see Popa |23| . de la Harpe and Jones ^Hl and 
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Munemasa and Watatani 20 . In the last two papers one construct complex Hadamard matrices 
not of Sylvester type when n is a prime number such that n = ±1 (mod 4). A little later Haagerup 
16 obtained the first example of a 6-dimensional matrix which is not a solution of the system of 
equations 

In this paper we make use of a few analytic techniques from the operator contraction theory 
and the factorization of unitary matrices to obtain a convenient reprezentation of unitary matri- 
ces of arbitrary order n that leads us easily to a system of (n — 2) 2 trigonometric (or equivalently 
polynomial) equations whose solutions give all the complex Hadamard matrices of order n. Our ap- 
proach is also useful for finding real Hadamard matrices, being complementary to the combinatorial 
approach almost exclusively used until now. 

The paper is organized as follows: in Section 2 the equivalence of the complex Hadamard 
matrices is reviewed. In Section 3 a theorem showing the existence of the complex Hadamard 
matrices for every integer n is stated and an upper bound on the number of continuum solutions 
is obtained. Section 4 contains an one-to-one parametrisation of unitary matrices written as block 
matrices and in the next Section an application of the obtained formulae is given. In Section 6 an 
other parameterisation of unitary matrices is given under the form of a product of n diagonal phase 
matrices interlaced with n — 1 orthogonal matrices each one generated by a real vector from R n . 
This form is convenient because it leads to a simpler form for the moduli equations and in the same 
time we consider it more appropriate for designing software packages for solving these equations. 
In Section 7 we show how to derive the moduli equations as trigonometric equations and give a few 
particular solutions for n = 6. In Section 8 the problem is reformulated as an algebraic geometry 
problem and we show that the parameterisation of Hadamard matrices can produce interesting 
examples for many problems currently under study in this field. The paper ends with Conclusions. 

2 Equivalence of complex Hadamard matrices 

Complex n-dimensional Hadamard matrices being unitary matrices whose entries have modulus 
the natural class of looking for complex Hadamard matrices is the unitary group U(n). 

The unitary group U(n) is the group of automorphisms of the Hilbert space (C n , (•, •)) where 
(•,•) denotes the Hermitian scalar product (x,y) = J2l=i ^iVi an d the bar denotes the complex 
conjugation. If A n £ U(n) by A* n we denote the adjoint matrix and unitarity implies A* n A n = 
A n A* n = I n . It follows that det A n = e tv , where <~p is a phase, and dimn U(n) = n 2 . 

Because in any group the product of two arbitrary elements is again an element of the group 
there is a freedom in choosing the "building" blocks to be used in a definite application. In the 
case of a complex Hadamard matrix the multiplication of a row and/or a column by an arbitrary 
phase factor does not change its properties and consequently we can remove the phases of a row 
and column taken arbitrarily. Taking into account that property we can write 

A n = d n A n d n —\ (2) 
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where A n is a matrix with all the elements of the first row and of the first column positive numbers 
and d n = (e l(fil , . . . , e llpn ) and d n -\ = (1, e llpn+1 , . . . , e 1 ^ 2 ™- 1 ) are two diagonal phase matrices. In 
the following we will consider that A n = A n , i.e. A n will be a matrix with positive entries in the 
first row and the first column. 

Since a unitary matrix is parameterised by n(n — l)/2 angles and n(n + l)/2 phases the 
above equivalence relation tell us that the number of remaining phases is n(n + l)/2 — (2n — 1) = 
(n — l)(n — 2)/2, and so the number of free real parameters entering a unitary matrix is reduced 
from n 2 to n 2 — (2n — 1) = (n — I) 2 . 

Secondly we can permute any rows and/or columns and get an equivalent unitary matrix. This 
procedure can be seen as a multiplication of A n at left and/or right by an arbitrary finite number of 
the simplest permutation unitary matrices Pij, i ^ j, i,j = l,...,n, whose all diagonal entries but 
an and ajj are equal to unity, an = ajj = 0, = dji = 1, i ^ j and all the other entries vanish. 
Both the diagonal phase and permutation matrices generate subgroups of the unitary U{n) group; 
so we may consider them as gauge subgroups, i.e. any element of U(n) is defined modulo the action 
of a finite number of the above transformation, which has as consequence a standard representation 
for unitary matrices. We consider that the group generated by the above two subgroups deseves to 
be independently studied since its orbit structure could shed light on many important issues from 
information theory and stochastic matrices. 

The above two equivalence conditions are those found by Sylvester [21] for the Hadamard 
matrices, but in fact they are valid for U(n) which is invariant with respect to the product of an 
arbitrary number of the above transformations. 

Besides for Hadamard matrices we will not distinguish between A n and its complex conjugated 
matrix A n , the complex conjugation being equivalent to the sign change of all phases (fi — > —ipi 
entering the parametrisation. More generally we shall consider equivalent two matrices whose 
phases can be obtained each other by an arbitrary non-singular linear transformation with constant 
rational coefficients. As we will see later the complex Hadamard matrices depend in general on a 
number of arbitrary phases and the above condition says that we will consider only the most general 
form of the solution and not those particular forms obtained by prescribing definite values to the 
(arbitrary) phases entering the parameterisation. In this sense we can say that there is only one 
complex Hadamard matrix of order 4, that found by Hadamard |17j . all the others, including those 
with all entries real numbers, being particular cases of the complex one. Other authors speak in 
this case of non-equivalent or a continuum of solutions ^S] • We consider that the above conditions 
are the only a priori equivalence criteria we can impose on Hadamard matrices, i.e. will consider 
equivalent any two matrices that can be made equal by applying them a finite number of the above 
transformations. 
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3 Existence of complex Hadamard matrices 



The parameterisation of a unitary matrix by the moduli of its entries is very appealing, and in the 
case of Hadamard matrices compulsory, although it is not a natural one in the general case. A 
natural parameterisation would be one whose parameters are free, i.e. there are no supplementary 
restrictions upon them to enforce unitarity. In this sense natural parameterizations are the Euler- 
type parameterisation by Murnagham L21 : , or that found by us [TT] . 

The problem we rose in was to what extent the knowledge of the moduli |ay| of an n x n 
unitary matrix A n = (a^) determines A n . Implicitly we supposed that A n is parameterized by n 2 
independent parameters. But from what we said before we know that we may ignore 2n — 1 phases 
entering the first row and the first column and consequently the number of independent parameters 
reduces to (n — l) 2 , that coincides with the number of independent moduli implied by unitarity. If 
we identify the parameters to the moduli they will be lying within the simple domain 

D = (0,1) x ... x (0,1) = (0,1) ( ™~ 1)2 

where the above notation means that the number of factors entering the topological product is 
(n — l) 2 . We excluded only the extremities of each interval, i.e. the points and 1 that is a 
zero measure set whitin U(n) and has no relevance to the parameterisation of complex Hadamard 
matrices. 

Thus, in principle, we can parameterise an n x n unitary rephasing invariant matrix by the 
upper left corner moduli; we exclude the moduli of the last row and of the last column since they 
follow from unitarity. Nothing remains but to check if the new parameterisation is one-to-one. A 
solution to the last problem is the following: start with a one-to-one parameterisation of U(n) and 
then change the coordinates taking as new coordinates the moduli of the (n — l) 2 upper left corner 
entries (and 2n— 1 ignorable phases). Afterwards use the implicit function theorem to find the points 
where the new parameterisation fails to be one-to-one. The corresponding variety upon which the 
application is not a bijective one is given by setting to zero the Jacobian of the transformation. One 
gets that generically for n > 4 the unitary group U(n) cannot be fully parametrised by the moduli 
of its entries, i.e. for a given set of moduli there could exist a continuum of solutions, but this 
negative result is good for the parameterisation of Hadamard matrices by decreasing the number 
of independent solutions taking into account the equivalence conditions discussed in the previous 
section. 

n(n — 3) 

If the moduli are outside of the above variety an upper bound for the multiplicity is 2 2 . 
However in the case of Hadamard matrices the equivalence constraints reduce this number to lower 
values than the above upper bound. The bound is saturated for n = 3 when there is essentially only 
one complex matrix, i.e. for given moduli values for the first row and column entries compatible 
with unitarity, the sole freedom is an arbitrary phase. If we denote the relevant squared moduli by 
mi , ni2 , 7773 , 777,4 and the phase by tp then the compatibility condition has the form 

— 1 < COS If = ( — 1 + 2777,1 — 777 2 + 7772 + ™3 + "74 — 777i7772 — 777i7773 — 77727773 — 
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2m\m\ — m\m2rn^m\m^) /2^J mirn2rn^{l — m\ — m^Hl ~~ m i ~ m 3) < 1 

This is also the necessary and sufficient condition which the squared moduli rrii, i = 1,...,4, 
have to satisfy in order to obtain a unistochastic matrix from a general doubly stochastic matrix. 
Because unitary matrices of arbitrary dimension do exist and on the other hand the number of 
independent essential parameters of a U(n) matrix is (n — l) 2 the following is true: 

Theorem 1 Suppose (xi, . . . , x n 2) is a co-ordinate system on the unitary group U(n) consisting 
of n(n — l)/2 angles each one taking values in [0, 7r/2] and n(n + l)/2 phases taking values in 
[0, 2ir). By discarding 2n — 1 non-essential phases the number of co-ordinates reduces to (n — l) 2 , 
(xi, . . . , £( ra _i)2), that coincides with the number of independent moduli (mi, ■ ■ ■ , ni(n-l) 2 ) implied by 
unitarity. Taking as new co-ordinates the moduli mj, i = 1, . . . , (n — l) 2 , the new parameterisation 
is generically not one-to-one for n > 4, the non-uniqueness variety being obtained by setting to zero 
the Jacobian of the transformation 

<9(mi,...,m (ra _ 1)2 ) _ q 
d(xi, . . . ,x (n _i )2 ) 

n(n — 3) 

Outside this variety the number of discrete solutions N s satisfies 1 < N s < 2 2 and on the 
variety described by J3]j there is a continuum of solutions. In the special case of complex Hadamard 
matrices all the solutions are given by the system of trigonometric equations 

mf(xi,...,xr n _ 1)2 ) = - , i = 1, . . . , (n - l) 2 (4) 

Suppose we know the irreducible components of the variety (pj) and let r(n) be the rank of the system 

in every irreducible component, then every solution of OJ) in such an irreducible component will 
depend upon (n— l) 2 — r(n) arbitrary parameters and the number of (continuum) solutions satisfies 

1 < N s < 2 r ( n )~ 1 ^ n ( n ^ 1 )/ 2 . 

Proof. In the general case Eqs.(@J) have the form 

m\(x\, . . . ,x 2 n _ 1) ) = a,i, where a% G (0, 1) , i = 1, . . . , (n - l) 2 (5) 

The parameters a, generate a doubly stochastic matrix. The Eqs.©, as we will see later, are 
trigonometric equations in our parameterisation, and consequently the multiplicity of the solutions 
may arise from the two possible phase solutions for all values of sine or cosine functions that 
satisfy (fKj). The number of independent phases is (n — l)(n — 2)/2 and taking into account that 
we consider A n and A n as equivalent matrices, condition which halves the number of solutions, the 
above bound for N s follows. A similar argument establishes the upper bound for the number of 
continuum solutions. 

For n = 3 the Jacobian is positive and 1 < N s < 1, which implies the existence of one complex 
matrix irrespective of the values aj, compatible with unitarity. 
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It is easily seen that the equations which correspond to the first row and the first column entries 
have a unique solution and the number of equations reduces to (n — 2) 2 . Indeed, because these 
entries are positive we can take the following parameterisation in terms of 2n — 3 angles, e.g. for 
the first row 

(an, . . .,a ln ) = (cosxi,sinxi cos\2, sinxn-i) 
and similarly for the first column. The Eqs.(J5J) give the unique solution 

cos Xk = — t—- i , k = 1, 2, . . . , n — 1 

iKa-*) 

where = | a-ifc | 2 , k = 1, 2, . . . , n — 1. In the case of Hadamard matrices one gets 

cog Xfc = -7 \ , , fc = 1,2, ... ,n — 1 
yn + 1 — k 

and the same solution for the angles parameterising the first column. In this way the number of 
equations reduces to (n — l) 2 — (2 n — 3) = (n — 2) 2 and the upper bound for the continuous solutions 
may be written as 1 < N s < 2 r ( n ) _1 ~( n ~ 2 )( n ~ 3 )/ 2 ; where r(n) is the rank of the reduced system. 
Even so the number of equations grows quadratically with n which shows that even for moderate 
values of n the problem is not easy to solve. 

In conclusion we have a system of trigonometric equations whose solutions will give all the com- 
plex Hadamard matrices, but to get effective we have to start with a one-to-one parameterisation 
of unitary matrices in order to find the explicit form of the (n — 2) 2 equations and try to solve 
them. In the following Section we will provide one of the two parameterisations of unitary matrices 
that we will use in the paper. 



4 Parameterisation of unitary matrices 

The aim of this section is to provide a one-to-one parameterisation of unitary matrices that will be 
useful in describing the complex Hadamard matrices. We shall present two such parameterisations 
and for the the first one we follow closely our paper showing here only the points which are 
important in the following. The algorithm we provide is a recursive one, allowing the parameter- 
isation of n x n unitary matrices through the parameterisation of lower dimensional ones. The 
parameterisation will be one-to-one and given in terms of a(n) angles taking values in [0,7r/2] and 
f{n) phases taking values in [0, 2tt) such that the application 

A n (A n G U(n),A n Al = I n ) - E = (0, vr/2) a W [0, 2tt)^ C R n * 

is bijective. Always in the following the ends of the interval [0, 7r/2] will be obtained by continuation 
in the relevant parameters, if necessary. 



8 



The starting point is the partitioning of the matrix A n 6 U (n) in blocks 

For definiteness we suppose the order of A is equal to m with m < n/2. The blocks entering © 
are contractions as follows from unitarity 

AA* + BB* = I m , A* A + C*C = I m , CC* + DD* = J n _ ro (7) 

where in the following Ij. denotes the x k unit matrix. Suppose we know the contraction A, then 
the problem reduces to finding the B, C and D blocks such that A n should be unitary. In other 
words the problem is: knowing a contraction A of side m how we can border it for getting a unitary 
nxn matrix A n ? For solving this problem we shall make use of the theory of contraction operators. 

An operator T applying the Hilbert space TL in the Hilbert space TL' is a contraction if for 
any v £ TC, \\Tv\\n' < |M|-h, i.e. ||T|| < 1, [22]. For any contraction we have T* T < In' and 
TT* < In and the defect operators 

D T = (I n - T* T) 1 / 2 , D T * = (I w - TT*) 1 ' 2 

are Hermitean operators in TL and TL' respectively. They have the property 

TD T = D T * T, T* D T * = D T T* (8) 

Here we consider only finite-dimensional contractions, i.e. T will have in general n\ rows and 
columns. 

The unitarity relations ® can be written as 

BB* = D%, C*C = D\ 

According to Douglas lemma ^1] there exist two contractions U and V such that 

B = D A »U, and C = D A V 

Since we are looking for a parameterisation of unitary matrices, U* and V are isometries, i.e. they 
satisfy the relations 

TIT I* — T V*V — T 

If n is even and m = n/2, then U and V are unitary operators. Thus B and C blocks are given by 
the defect operators Da* , Da and two arbitrary isometries whose dimensions are m x (n — m) and 
(n — m) x m respectively. The last block of A n is given by the lemma 
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Lemma 1 The formula 

D = -VA*U + D v ,KDu 



establishes a one-to-one correspondence between all the bounded operators D such that 

A n - 



A D A *U 
VD A D 




is a contraction and all the bounded contractions K. 

See Arsene and Gheondea [2| for a proof of the general result when U, V and K are contractions, 
and further details. In our case U and V being isometries D is given by 

D = —VA*U + XMY (9) 

where X and Y are those unitary matrices that diagonalise the Hermitean defect operators Dy* 
and Djj respectively, i.e. 

X*D V *X = P, Y*DuY = P 

P is the projection 

P = 

and the matrix M entering @ has the form 

M = 

where A n _2m denotes an arbitrary (n — 2m) x (n — 2m) unitary matrix. See jllj for details. In 
the above formulae we supposed that the eigenvectors of the Djj and Dy* operators entering the 
matrices X and Y are ordered in the increasing order of the eigenvalues. 

Therefore the parameterisation of an n x n unitary matrix is equivalent to the parameterisations 
of four matrix blocks with lower dimensions than those of the original matrix, and consequently 
our task is considerably simplified. On the other hand the formulae ® and subsequent show that 
this procedure is recursive allowing the parameterisation of any finite dimensional unitary matrix 
starting with the parameterisation of one- or two-dimensional unitary matrices. Moreover the 
parameterisation of A n requires the parameterisation of an m x m contraction, of two isometries 
U and V and of an (n — 2m) x (n — 2m) unitary matrix. In our papers [111 I12j we considered only 
the case m = 1 as the simplest one, however the case m > 1 may be useful in the study of complex 
Hadamard matrices. 

For what follows we treat again the case m = 1, i.e. A is the simplest contraction, a com- 
plex number whose modulus is less than one, because we found the form of the matrices X and 
Y for arbitrary n. Since V is a (n — l)-dimensional vector the isometry property allows us to 
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parametrise it as V = (cos xi, sin xi cos X2, ■ ■ ■ , sin Xi ■ ■ ■ si n Xn-2) where t denotes transpose. V 
is the eigenvector of Dy* corresponding to the zero eigenvalue. Indeed from the relations (JHJ) we 
have 

D v *V = VD v = 

showing that V is the eigenvector of Dy* corresponding to the zero eigenvalue. Thus the problem is: 
how to complete an orthogonal matrix X knowing its first column (row) such that no suplementary 
parameters enter. The other columns of this matrix we are looking for will be given by the other 
eigenvectors of Dy*. One easily verifies that Dy* is a projection operator such that the other 
eigenvalues equal unity. Indeed the folowing holds 



Lemma 2 The orthonormalised eigenvectors of the eigenvalue problem 

D v *v k = X k v k , k = l,...,n-l 

are the columns of the orthogonal matrix X £ SO(n — 1) and are generated by the vector V as 

( cos xi \ 

sin xi cos X2 

Vl 



\ sinxi--- sinXn-2 J 



and 



d 



Vk+l 



7T 



Xk-i 



), k = 1, . . . , n 



dxk 2 
where in the above formula one calculates first the derivative and afterwards the restriction to ir/2. 

In a similar way one finds Y; see |13j for a proof. 

In the case of n x n Hadamard matrices whose elements of the first row and of the first column 
are positive numbers a±j = dji = j = 1, . . . , n, X has the form 



fn-X 
1 



n-2 

n-l 



X 



n-l y/(n-l)(n-2) 
1 1 



n— 3 
n-2 








y/(n-X)(n-2) v / (n-2)(n-3) 



n— 4 
n— 3 



n-l y/(n-X)(n-2) 
1 1 



y/(n-2)(n-3) 



y/(n-3)(n-4) 



\ v^T y/(n-X)(n-2) y/(n-2)(n-3) y/(n-$)(n-t) 
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and Y = X , where t denotes the transposed matrix. 

In this way all the quantities entering formula © are known and the parameterisation of A n 
can be obtained recursively starting with the known parameterisation of 2 x 2 unitary matrices. 

When the block A is one-dimensional, i.e. a simple number equal to the term V A* U 

entering Eq.® has the form J where J is the (n — 1) x (n — 1) matrix whose each of entries 

is +1, which appears in many constructions of real Hadamard matrices; see Agaian pQ. 



5 Application 

In the following we will use Eq. Q to generalize to the case of complex Hadamard matrices the tries 
used by Sylvester [21] and Hadamard JJ] for constructing complex Hadamard matrices. We take 
n an even number, n = 2 m, and we suppose that we know a parameterisation of the A block which 
is unitary and whose order is m. In that case B and C blocks are also unitary matrices of order m 
and we consider them normalized as A A* = B B* = C C* = I m . From © we have D = —C A* B 
and the following matrix 

J_(A B \ 

V2\C -CA*B J 

will be unitary by construction. In general the above matrix will not be Hadamard even when 
A, B and C are, as the simplest example shows; this happens only when either C = A or B = A. 
Since the second case is obtained by transposing the matrix of the first one, as long as B and C 
are arbitrary, we will consider only the matrix 

° X (10) 

which is the elementary two-dimensional array that will be used in the construction of more com- 
plicated arrays of Hadamard matrices. In the following we suppose that A and B are complex 
Hadamard matrices of size m each one depending on p > 0, respectively, q > free phases, i.e. 
(110(1 is a complex Hadamard matrix of size 2 m. Now we make use of Hadamard's trick to get a 
Hadamard matrix depending on p + q + m — 1 arbitrary phases. Indeed we can multiply B at left by 
the diagonal matrix d = (1, e l . . . , e l ^ m_1 ) without modifying the Hadamard property. In this 
way Hadamard obtained a continuum of solutions for the case n = 4. We denote B\ = d ■ B and 
then the matrix 




J_ / A B l 
V2\ A -Br 



11) 



will be unitary and Hadamard depending on p + q + m — 1 parameters. From (|llj) we obtain in 
general two non-equivalent 2 m x 2 m Hadamard matrices when B 7^ B* . In this case Eq. (jlljl is a 
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realization and the second one is given by B\ —* B 2 = d ■ B* . The above procedure can be iterated 
by taking the matrix ()1U|) as a new A block obtaining a Hadamard matrix of the form 



2 



1 



(A B C D\ 

A -B C -D 

A B -C -D 

\ A -B -C D ) 



(12) 



which is a 4 m-dimensional array similar to Williamson array [2HJ ; and so on. In contradistinction 
to the Williamson array the A, B, C, D blocks satisfy no supplementary conditions, excepting their 
unitarity. Thus the following holds 

Proposition 1 If the m x m complex Hadamard matrices A, B, C, D depend on p, q, r, s arbitrary 
phases then there exists a complex Hadamard matrix of the form which depends on p + q + r + 
s + 3(m — 1) arbitrary phases. 

We notice that the elementary array IjlOjl is different from the Goethals-Seidel one |15j that 
appears in the construction of real Hadamard matrices and which has the form 



being the relation AB* = BA*. We consider that the form ijTUjl could also be useful for the study 
of orthogonal designs and real Hadamard matrices it being in some sense complementary to the 
above form. 

As an application of the formula (|12|) we consider the following case: an = 012 = 021 = —a 22 = 
hi = &12 = en = C12 = d u = d 12 = l/y/2 and b 21 = -b 22 = e ts /V2, c 2 i = -c 22 = e lt /V2, 
d 2 \ = —d 22 = e m /\/2 where the notation is self-explanatory, and we obtain an eight-dimensional 
Hadamard matrix depending on three arbitrary phases s, t, u. 

When A = B, Eq. ljlfljl can be written as 



where e = — 1, i.e. the first factor is the Sylvester Vandermonde matrix of the second roots of 
unity, and (8> is the ordinary Kronecker product, A (g> B = [aijB]; of course the first factor can be 
any complex Hadamard matrix of order m. Now we want to define a new product the aim being a 
more general construction of Hadamard matrices. Let M and N be two matrices of the same order 
m whose elements are matrices Mij of order n and respectively N^i of order p. The new product 
denoted by (g> is given as 




The above array is not unitary even when A and B are, the suplementary condition for unitarity 




(13) 



Q = M®N 
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which is a matrix of order mnp, where 



k=m 

Qij = M ik ® N kj 

k=l 

We will use here the above formula only in the case: M = rriij where rriij are complex scalars, 
not matrices and N is an arbitrary diagonal matrix N = (iVn, • • • , N mm ) where Na ar matrices of 
order p obtaining 



Q 



I miiiVn 



rn 



N n 



(14) 



nun 1 ' inrn 



Thus the following is true. 

Proposition 2 If the matrices M and Na, i = 1, . . . ,m, are Hadamard so will be the matrix Q 
given by Eq.^T^j. 

The order of Q is mp and the formula (|14|) is new even for real Hadamard matrices. This form is 
the most general array we have obtained and in some sense ()14j) is the natural generalization of 
Williamson arrays to the case of complex Hadamard matrices. 

If in the above relation we take m\\ = m\2 = 77121 = — "2-22 = l/v2 and N\\ = A and ./V22 = B, 
then Eq.JHJ) reduces to Eq. (fTn|) . 

Example 1 // now rriij are the same as above and 







( 1 1 


1 












1 1 


-1 


-1 








Nu = - 
2 


1 -1 


_ e is 












V 1 - 1 










complex four- dimensional Hadamard matrix and 










( 1 ° 


\ 


/I 


1 


1 


1\ 




e u 





1 


1 


-1 


-1 


N 22 = \ 


e 


in q 


1 


-1 -e^ 






^ 


e iv ) 


V 1 


-1 e*» 





we obtain an eight- dimensional matrix depending now on five arbitrary phases s,t,u,v,y instead of 
three as in the preceding example obtained by using the Williamson-type array (12). 



Thus the following holds. 
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Proposition 3 If M,Ni, i = l,...,m are m x m and respectively, n x n- dimensional complex 
Hadamard matrices depending on m, respectively, rii, arbitrary phases, then there is an array of 
the form (14) that depends on 



m + n\ + (m — 1) ^ 

i=2 



free phases. 



The above example shows the necessity for getting upper and lower bounds on the number 
of arbitrary phases entering a Hadamard matrix of size N. Taking into account the standard 
decomposition of any integer under the form N = p qi . . . p^ 1 , where p\ < . . . < p m are primes and 
qi . . . q m their respective powers, we may use the above Proposition 3 for obtaining lower bounds 
on the number of free phases, that we shall denote it by (f{N). Since until now does not exist 
an example of a Hadamard matrix of size N with N prime which depends on free phases, in the 
following we will consider the normalization tp(N) = 0, for N prime. Thus the following holds. 

Theorem 2 Let N = p^ 1 be the power of a prime p\, with q > 2. Then a lower bound for ip(pl 1 ), 
the number of free phases entering the parameterization of the N x N complex Hadamard matrix, 
is given by 

^(p? i ) = i+[(pi-i)(?i-i)-i] P r 1 

IfN = pf...p%?= pfNi then <p(pfN{) is given by 

ifipfm) = 1 + [(pi - l^iiVi - Pl ]pf - l + <p{N{)p? 
Proof. Making use of Proposition 3 we find the recurrence relation 

M 1 )=pMpT~ 1 ) + (pi-i)(pT' 1 -i) 

with the initial condition <p(pi) = and the solution follows. 
In the second case the recurrence relation writes 

Mm) = pMpr'Ni) + (pi - w'jvx - 1) 

and the initial condition can be taken as 

¥>(piJVi) = pMN x ) + (p! - l)(iVi - 1) 

and the solution follows. The above recurrence relation allows us to obtain lower bounds for any 
integer N under the form 

ip{pf ...p q ™) = l + [(pi - l)q lP f ...p q ™- Pi]p\ L ~ 1 + 

pf {1 + [(p 2 - \)q 2P f . .. P <% - P2 ]p q r 1 + 
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pf{l + [(p 3 - l)q 3 pf . ..p&* - P3 ]pf~ 1 + 
pf{l + ... +P^T 1 1 {1 + [(Pra ~ l)q m - PmK m " 1 } + 

pS{1 + liPrn - l)(q m " 1) " Mp^ 1 } ■ • •} 
We give now a few examples. 

Example 2 If N = p^p^ then the lower bound for (fiPiP^), the number of free phases entering 
the parameterization of the N x N complex Hadamard matrix, is given by 

<p(pVp¥) = i + (pi - i)<rf ~ x pf + [(p 2 - 1)(92 - 1) - iK^r 1 (is) 

Numerical examples: (y9(2 3 ) = 5, v?(2 4 ) = 17, <p(6) = 2, (^(3 2 ) = 4,^(2 2 3 2 ) = 49, etc. 



6 An other parameterisation of unitary matrices 



In the following we will shortly present another parameterisation of unitary matrices |13j under 
the form of a product of n diagonal matrices containing phases interlaced with n — 1 orthogonal 
matrices each one generated by a real vector v G R™. This new form will be more appropriate 
for design and implementation of the software packages necessary for solving the equations Q for 
arbitrary n. 

We have seen in Section 2 that we can write any unitary matrix as a product of two diagonal 



matrices of the for d n = (e 



,e lLPn ) with ifj £ [0,27r), j = 1, ...,n arbitrary phases and a 



unitary matrix with positive elements in the first row and the first column. We make also the 



notation d\ 



(In— k 5 ^ 



ii/'i 



k < n, where l n _fc means that the first (n — k) diagonal 



entries equal unity, i.e. it can be obtained from d n by making the first n — k phases equal zero. 
These diagonal phase matrices are the first building blocks in our construction. Other building 
blocks that will appear in factorization of unitary matrices A n are the two-dimensional rotations 
which operate in the i,i + 1-plane of the form 



J, 





V o 







cos ( 
sin i 



— sm Vi 
cos Qi 







(16) 



In-i-1 / 



The factorization idea comes from the well known fact that U{n) acts transitivly on the n- 
dimensional complex sphere S2n-i £ C n , and explicitely from the coset relation 



S2n-i = coset space U(n)/U(n — 1) 
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A direct consequence of the last relation is that we expect that any element of U (n) should be 
uniquely specified by a pair of a vector v G S271-1 an d of an arbitrary element of U(n— 1). Thus 
we are looking for a factorization of an arbitrary element A n G U(n) in the form 



A 



n 




where B n G U (n) is a unitary matrix whose first column is uniquely defined by a vector v G S2n-i> 
but otherwise arbitrary, and ^n-i is an arbitrary element of U(n — 1). Iterating the previous 
equation we arrive at the conclusion that an element of U (n) can be written as a product of n 
unitary matrices 

A n = B n ■ . . . -B" -1 



where 

^rt— k 




Bfr, k = 1, . . . ,n — 1, are k x k unitary matrices whose first columns are generated by vectors 
bk G S2fc-i; for example is the diagonal matrix (1, . . . , 1, e liPn ( n + 1 '>). 

The still arbitrary columns of B^ will be chosen in such a way that we should obtain a sim- 
ple form for the matrices B^~ k , and we require that Bk should be completely specified by the 
parameters entering the vector bk and nothing else. 

If we take into account the equivalence considerations of the Section 2 then B n (B n -k) can De 
written as 

Bn = d n B n 

where the first column of B n has non- negative entries. 

Denoting this column by v\ we will use the parameterization 



v\ = (cos 0\, cos 02 sin 0\,..., sin 0\ . . . sin n -i) t 



where 0{ G [0, vr/2], i = 1, . . . , n — 1. Thus B n will be parameterized by n phases and n — 1 angles. 
According to the above factorization B n is nothing else than the orthogonal matrix generated by 
the vector v\ and its form is given by Lemma 2 with n — > n + 1. Thus without loss of generality 
B n = d n O n with O n G SO(n). In this way the factorization of A n will be 

A n = d n On 4-1 <?n-l ■ ■ ■ d^O^d^ 1 I n (17) 

where 0\_ k has the same structure as B^_ k , i.e 



O k 



n—k 



o 



n—k 



and = (l,...,l,e^,...,,e^ 



k 



The orthogonal matrices O n can be factored in terms of Jj,j+i as follows. 



17 



Lemma 3 The orthogonal matrices O n ( 0^_ k ) at their turn can be factored into a product of 
n — 1 (n-k-1) matrices of the form J%i+\; e.g. we have 

@n — Jn—X,n Jn—2, n-l ■ ■ ■ Jl,2 
where Ja+i are nx n rotations introduced by Eg. Ub\) . 

In this way the parameterisation of unitary matrices reduces to a product of simpler matrices: 
diagonal phase matrices and two-dimensional rotation matrices. For more details see our paper 
J13J. Now we propose a disentanglement of the angles and phases entering each "generation" and 
denote the angles by latin letters, e.g. those that parameterize O n will be denoted by a\, . . . , a n _i, 
the angles that parameterize O n _ 1 , by b\, ... , 6 n -2, etc., the last angle entering 0^~ x by z\. The 
phases will be denoted by Greek letters; e.g. the phases entering d\ will be denoted by ati, ■ ■ ■ , a n , 
those entering d n _i by (3%, . . . ,/? n -i> etc. The above factorization will be used in the next section 
for obtaining the equations for the moduli of the matrix elements. 



7 Explicit equations of the moduli 

Our choise for the orthogonal vectors in Lemma 2 was such that the resulting matrix should have as 
many zero entries as possible. Thus O n has (n — l)(n — 2)/2 zeros in the right upper corner and the 
entries of the Hadarmard matrix will get more and more complicated when going from left to right 
and from top to bottom. We will start using the form (|17|) of the unitary matrix and then d n = I n . 
Since the first column has the form an = 1/ 1 \pn, i = 1, . . . , n and dj l _ 1 = (1, e ia , e iai , . . . , e tan - 2 ) 
the product O n d\_ x is 



( ± 

1 

\fn 
1 


In— 1 „ia 

-y — e 

i/n(n-l) 

£ ia 




V n-l e 





In— 3 „ia 2 








\ 





\fn 


yjn(n-l) 


y/(n-l)(n-2) 


V «-2 e 


1 


e ia 




e ia 2 


e ia„_3 


— e ia n-2 


\fn 


i/n(n-l) 

e ia 


7(n-l)(n-2) 
e ia l 


V(n-2)(n-3) 

e™2 




V2 

e ia n-2 




y/n(n-l) 


y/(n-l)(n-2) 


V(n-2)(n-3) 


V6 


V2 J 



where a, «j, i = 1, . . . , n — 2 are n — 1 arbitrary phases. 
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The next building block 0\_ x A 2 n _ 2 will have the form 



/ 1 









cos a 
sin a cos a\ 



sin a e 



cos a cos a\ e 



%p 



\ sin a... sin a n _3 cos a sin a\ ■ ■ ■ sin a n _3 e 



»0 



cos a n _3 e 



iPn-3 



(19) 



in terms of n — 2 phases f3,f3x,..., /3 n -3 and n — 2 angles a, Oi, . . . , a n _3, and so on. 

It is easy to see that the first two columns of the product of matrices (18) and (19) does not 
change when multiplied by C^_ 2 d^_ 3 ; however the first row does. If the angles entering 0\_ 2 are 
denoted by b, bi, . . . , 6 n _4 and the phases are 7, 71, ... , 7 n -4, etc., then the entries of the first row 
are 



012 



n 



1 



■ cos a e 



n 



«13 



11 



1 



■ sin a cos 



n 



6e i(a+/3) ,.... 



flln-l 



-1 



,n-l 



• sin a sin b . . . cos z e 



i(a+/3+...w) 



where z and w are the last angle and phase respectively. Since we use the standard form of 
Hadamard matrices, i.e. the entries of the first row and of the first column are positive and 
equal l/\/n, the above equations imply 



a = f3 



lo = 7r; cos a 



n 



1 



cos b 



, cos z 



n 



1 

71 



We substitute the above values in Eq, ()17|) and find a complex n x n matrix depending on (n — 
l)(n-2)/2 phases ot\, . . . , a n - 2 , fh, • • • , "4>l and (n-2)(n-3)/2 angles ai, . . . , a n _ 3 , 61, . . . , y\, i.e. 
(n — 2) 2 parameters which have to be found by solving the corresponding equations given by the 
moduli. The first simplest entries of the unitary matrix have the form 



1 



n 



a-22 



(n - l)y/n 



n 



■ cos ai e 



UXl 



Ofc2 = — 



(n - l)y/n 



+ 



' n 



cos a\ e 



n 



v/(n - l)(n 



+ ...+ 



' n — k 
n — k + 1 

1 



sin ai . . . sin cifc_2 cos a^-i e 



(n - l)-v/n 



+ 



cos ai e 



(.HI 



sm ai . . . cos afc_2 e 
^(n-/c + 2)(n- A; + l) 

fc = 3, . . . , n — 1 

sin ai cos 61 e 1 ^ 1 " 1 "^ 1 ^ 



(20) 



n 



1 U(n-l)(n-2) v /( n _ 2 )(n-3) 



+ ... 
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+ (-l) fe - 1 J n - k - sin ai sin h . . . cos 1(h)! e *(°i+ft+-+*(*)i) j , etc. 
y n — A; + 1 / 

where and X(k) denote the letters for angle and respectively phase corresponding to index k 
and the signs in the last bracket alternate. 

The matrix elements get more complicated when going from the upper left corner to right 
bottom corner. The entries 022,032 and 023 lead, for example, to the following moduli equations 

2 2 

(n — 2) cos ai H -= cos a\ cos a.\ — 1 = 

Jn 



sin ai ((n — 3) sin a\ cos 2 02+ 

ruh u > i — =— - cosai cos(ai — 02) I — sin ai 

n - 1 V V n / / 

sin ai ^(n — 3) sin a\ cos 2 61+ 



In — 3 / cos a2 

■ cos 02 — -j= cos a\ cos(ai — ai) - sin ai =0 (21) 



In -3 ( cos(ai+/?i) 

■ cos 61 t= h cos ai cos ft - sin ai =0 



y n — 1 V y^n 

and so on. The form of the last two equations was obtained after the elimination of the term contain- 
ing cosai cosai by using the first equation (21), i.e. we work in the ideal generated by the moduli 
equations. It is easily seen that the other equations contain as factors sin 02, • • • , sina n _2, sin 61, . . . , etc.. 
Thus a particular solution can be obtained when 

sin a\ = 

which implies a\ = 0, it, and from the first equation (21) we get 

An- 3)Vn 

cosai = ± 

2 

It is easily seen that the above equation has solution only for n = 2,3,4; for n > 5 the factor 
sinai will be omitted from Eqs.(21) because then a\ 7^ 0,7r. When n = 2 we obtain a\ = vr/4, so 
022 = —l/y/2. If n = 3, then a\ = 3ir/2 and from the first Eq.(20) one gets 

1 i 1 27Ti 

a22 = -27i + 2 = 7r 3 ' etc - 

The case n = 4 leads to a\ = ir which gives 

X e i(oL2+0i) 
a 2 2 = -023 = -032 = ^ aIld °33 = -«34 = ^ 
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After the substitution oli + (3\ = t one finds the standard complex form of the 4x4 matrix found 
by Hadamard. To view what is the origin of the phase a<i + (3\ we have to look at the moduli 
equations. They have the form 

2 cos a\ + cos a\ cos a% — 1 = 

sinai(cosa2 — 2cosai cos(ai — 02)) = 

sin ai (2 cos a\ cos f3\ — cos(ai + j3\)) = 

cos 2ai cos(«i — 02) cos j3\ + cos a\ cos(«2 + fix) + sin(ai — 02) sin (3\ = 

and we see that the above system splits into two cases. In the first case, when sinai = 0, the rank 
of the system is two which explains the above dependence of 033 on two phases and in the second 
case when sin a\ 7^ the rank is three and the dependence is only on one arbitrary phase. However 
in this case there is no final difference between the two cases. The solution of the above system 
is obtained directly but for n > 5 the problem is difficult and needs more powerful techniques. 
Particular solutions can be obtained rather easily e.g for n = 6 there is a matrix that has the 
property a^- = a^. 



There exists even a Hermitian matrix S 





( 1 


1 


1 


1 


1 


1 


\ 




1 


-1 


-1 


1 


i 


— i 




1 


1 


-1 


— i 


-1 


1 


i 






1 


1 


-1 


—i 


-1 


i 






1 


i 


1 


-1 


-1 


—i 






V 1 


—i 


i 


i 


—i 


-1 


/ 


matrix S 


= S* 














/ 1 


1 


1 


1 


1 


1 


\ 




1 


-1 


i 


/' 


—i 


—i 




1 


1 


—i 


-1 


1 


-1 


i 






1 


—i 


1 


-1 


i 


-1 






1 


i 


-1 


—i 


1 


-1 






V 1 


i 


—% 


-1 


-1 


1 


/ 



and so on. As we said before getting the most general form of a solution is not a simple task; for 
n = 6 we have 16 complicated trigonometric equations and we remind that the simpler (^Q) system 
was solved only for n < 8 equations. Thus new approaches are necessary and in the next Section 
we suggest such an approach: using methods from algebraic geometry. 
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8 Connection with algebraic geometry 

The Eqs.(21) can be transformed into polynomial equations by the known procedure 



sin a 



2x 

TT^2 ' 



cos a 



1-x 2 
1 + x 2 



such that we get from (21) 
Pi 



2 2 
(n - 3 + ~j=)x\ - 2(n - \)x\ + (n - 3 - -=] 



yj + (n - 3 - -j=)xj- 



2(n- l)x? + (n - 3 + -=) 



P2 



-(1 



d + ^)}^I + { 



n vn 



(1 



^2 ~~ Xl X 2 + (1 



n 



)x? + C\x\— 



L )x? + Cixi - (1 + -=) 
n Jn 



x\ - C 2 x\x\ - (1 



1 



n 



)xt+ 



d Xl + (l + -=)\yl + 



(1 + -L)xl + Cixi - (1 - 4=) 
n Jn 



1 



(1 + + Ci xi + (1 - -=) \ yi - 4(1 - xi)(l - x 4 2 ) ym + 



y/n 



y/n 



x\ - C 2 x 1 x 2 l - 



Ci xi + (1 - 



x| - C 2 xi x 2 , + (1 + — )x? + Cixi - (1 - — ) 
y/n Jn 



P3 



(1 



1 N x? + dxi + (i + -^; 



x 3 - C 2 xix| + (1 



1 



n 



)x\ + C\x\- 



(i - -\=)x\ + Cixi - (i + 



1 



x 3 - C 2 xi x 3 - (1 ^) x i + 



Cisi + (1 + -=)U? + 



.(1 * ^ + ^ + (1 * ) 



(l + 4=)aJi + Ci«i-(l-4=; 



where 



Cixi-(1- 



Ci 



1 



y3 2 -4(l + x?)(l-x|) yi y2 + 
1 



x 3 - C2 Xl x 3 + 
1 s 



x| - C 2 xi x| - (1 + — )x? + Cixi + (1 - _ 
\/n Jn 



y/n 

(n - l)(n - 4) 
y/(n- l)(n -IjJ' 



C 2 



2(n - l)(n - 2) 
y/(n-l)\n-^ 
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and the angles by the above transformation go to xi, X2, X3, . . . and the phases to yi, y2, 2/3, . . . 

From the matrices such as (18) one sees that the full set of the (n — 2) 2 equations contains 
square roots of almost all prime numbers < n so that not all the coefficients are rational and we 
have to look for solutions in a field Q(Vd) for some dGN. 

The polynomial equation p\ = defines an algebraic curve; however the most studied are the 
elliptic and hyperelliptic curves, i.e. those defined by an equation of the form y 2 = f p (x) where 
f p (x) is a polynomial of degree p. 

From pi = we get 

2 = (n-3-^)aj-2(n-l)sj + (n-3 + ^) = P l{xi ) 
Vl ~ (n-3+^ r )xf-2(n-l)x? + (n-3-^ r ) " P 2 (*i) 

which defines a meromorphic function. Its zeros and poles are 




'n + y/n - 2 



n — \ n 



and 

±1 



yn + 1 I n -y/n -2 

y/n-V \] n + y/n-2 



respectively that are simple, and the poles and the zeros are interlaced. Thus apparently the above 
equation is not hyperelliptic, however by the birational transformation 

Yi 



P 2 (xi) 
we get the equation 

Y? = -P 1 (x 1 )P 2 (x 1 ) 

which shows that the above curve has genus g = 3. For n > 5 the curve has no branch going to 
infinity since the highest power coefficient is negative and consequently the curve is made of three 
ovals. 

The polynomials p\ = p 2 = define a surface, p\ = p 2 = P3 = define a 3-fold, and so on. 
We consider that the study of these multi-fold varieties will be very interesting from the algebraic 
geometry point of view and their parameterizations could reveal unknown properties that may 
lead to a better understanding of the rational varieties. As we saw in Section 5 one can easily 
construct parameterizations of Hadamard matrices depending on a number of free phases at least 
for a non-prime n. That means that the set of the moduli equations has to be split in some sub-sets 
and for each such sub-set the solutions are in 5 1 (g> . . . (g) S\ , where k is the number of arbitrary 

k factors 

phases parameterizing the considered sub-set. But this could be equivalent to the existence of a 
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rational parameterization for the equations defining this sub-set. Unfortunately the best studied 
case and the best results are for algebraic curves; see |19| . Theorem 14, for a flavour of recent 
results. The study of surfaces, three-fold, etc. is at the beginning and until now the theory was 
developed only for the simplest varieties, the so called rationally connected varieties |19j . From 
what we said before one may conclude that the parameterization of complex Hadamard matrices 
could be an interesting example of the parameterization of meromorphic varieties, which could 
be a mixing between a rational parameterisation and a parameterisation of hyperelliptic curves. 
Thus the theoretical instrument for the parameterization of complex Hadamard matrices seems to 
exist, the challenging problem being the transformation of the existing theorems into a symbolic 
manipulation software program able to find after a reasonable computer time explicit solutions at 
least for moderate values of n. 



9 Conclusion 



All the results obtained for the complex Hadamard matrices can be used for the construction of 
real Hadamard matrices the only supplementary constraint being the natural one n = 4 m. We 
believe that the Hadamard conjecture can be solved in our formalism since unlike the classical 
combinatorial approach we have also at our disposal (n — l)(n — 2)/2 phases, and the problem is 
to guess the pattern of and tt taken by them. 

Conversely many constructions from the theory of real Hadamard matrices can be extended to 
the complex case. For example a complex conference matrix will be a matrix with an = 0, % = 
1, . . . , n and |ajj| = l/y^ such that 

WW* 

n 

It is not difficult to construct complex conference matrices, in fact it is a simpler problem than the 
construction of complex Hadamard matrices because the equations an = 0, i = 2, . . . , n — 1 imply 
the determination of 2(n — 2) parameters which simplify the other equations. 
We give a few examples: 
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where the second depends on two arbitrary phases. They are useful because if W n is a complex 
conference matrix then 





n 



In 



is a complex Hadamard matrix of order In. 

In this paper we have used convenient parameterisations of unitary matrices that allowed us 
getting a set of (n — 2) 2 polynomial equations whose solutions will give all the posible parameteri- 
sations for Hadamard matrices. Unfortunately the system is very complicated and only particular 
solutions have been found; thus from a pragmatical point of view the most important issue would 
be the design of software packages for solving these equations. 
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